Electric and magnetic dipole coupling in near-infrared split ring metamaterial arrays 
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We present experimental observations of strong electric and magnetic interactions between split ring res- 
onators (SRRs) in metamaterials. We fabricated near-infrared (1.4 /im) planar metamaterials with different 
inter-SRR spacings along different directions. Our transmission measurements show blueshifts and redshifts of 
the magnetic resonance, depending on SRR orientation relative to the lattice. The shifts agree well with a simple 
model with simultaneous magnetic and electric near-field dipole coupling. We also find large broadening of the 
resonance, accompanied by a decrease in effective cross section per SRR with increasing density. These effects 
result from superradiant scattering. Our data shed new light on Lorentz-Lorenz approaches to metamaterials. 

PACS numbers: 42.70.-a, 42.25.-p, 78.20.Ci 



Since the seminal work of Veselago and Pendry [1], many 
experimentalists have started to pursue optical materials with 
negative permittivity e and permeability n |2|]. The key moti- 
vation is the prospect of 'transformation optics', which allows 
arbitrary bending of electromagnetic fields, provided one has 
full control over e and fi. Particularly exciting examples are 
perfect lenses, that allow perfect sub-diffraction focusing JH, 
and 'cloaks' in which light passes an object without scatter- 
ing JH. Full control over e and /x requires 'metamaterials' of 
artificial nano-scatterers with electric and magnetic response, 
arranged in sub-wavelength arrays. The archetypical building 
block is the split ring resonator (SRR) consisting of a single 
cut metal loop with an inductive response. In recent years 
the field of metamaterials has made tremendous progress in 
shifting the resonant response from microwave to optical fre- 
quencies 0. |H H B Si • An important conceptual question is 
whether the effective response captured by e and /i is influ- 
enced by coupling between constituents. Coupling between 
SRRs in vertical ID stacks H9L 1 1 Qfl has attracted great atten- 
tion lately outside the scope of metamaterialSje.g., for mag- 
netic waveguides JHQj, 12, 131. antennas ffl . metamate- 
rial lasers B15I1 . and stereomaterials llal . Although constituent 
coupling might be anticipated to affect effective medium pa- 
rameters 11711 . measured effective responses have been at- 
tributed to single constituents in all experiments on metama- 
terial arrays to date. 

In this Letter we present the first measurements of strong 
constituent coupling in planar SRR metamaterial arrays. We 
fabricated and characterized SRR lattices with a magnetic re- 
sponse at A = 1.4 /zm JUla] in which we vary the spacing be- 
tween SRRs along different lattice directions independently. 
We observe large redshifts and blueshifts in the transmission 
resonances depending on SRR orientation relative to the lat- 
tices. We establish that in-plane electric-electric dipole cou- 
pling and out-of-plane magnetic-magnetic dipole coupling are 
strong competing interactions. We explain the shifts by a qua- 
sistatic electric and magnetic dipole coupling model 11 Oil , that 
enables us to determine the magnetic and electric polarizabil- 
ity of SRRs. Finally, we discuss the role of dynamic effects 
on the metamaterial resonance, which are evident in density- 
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FIG. 1: (Color) (a) We fabricated arrays of Au SRRs on glass with 
periodicities d x , y = 300 nm (SEM micrograph) and larger (inset). 
Each SRR has I = 200 nm, w = 80 nm, and SRR height 30 nm 
(±5 nm). (b) Transmission spectra for square SRR arrays with split 
width d = 80 nm (polarization along x). The magnetic resonance at 
1.4 fim blue-shifts and broadens with increasing density. 



dependent broadening and a saturation of the transmission. 

We have fabricated Au SRRs on glass substrates by elec- 
tron beam lithography and lift-off using PMMA resist 11311 . 
without any adhesive layers. We took great care to produce 
SRRs of identical dimension in arrays of different densities, 
using image analysis of SEM micrographs (see Fig. to 
overcome proximity effects. Based on |5||, our SRRs (200 nm 
base) are expected to have an LC resonance at 1.4 fim. Al- 
though driven by the electric field [18], we refer to the reso- 
nance as 'magnetic', consistent with literature [2]. To resolve 
the coupling strength between SRRs along the x (along SRR 
base) and y (along SRR arms) directions separately, we varied 
the pitches d x and d y independently between 300 nm and 550 
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FIG. 2: (Color)(a) Frequency of the magnetic resonance versus lat- 
tice spacing. The frequency blue-shifts when decreasing d y whether 
d y = d x (black squares) or not (blue circles, d x — 500 nm), while it 
red-shifts when decreasing d x (red triangles, d y = 500 nm). The 
inset shows raw spectra for d x — d y — 300 nm (black curve), 
d x = 500 nm, d y — 300 nm (blue curve), d x = 300 nm, 
d y — 500 nm (red curve), (b) FWHM of the magnetic resonance 
versus lattice spacing (color coding as in (a)). Curves are theory 
(electrostatic in (a), electrodynamic in (b)). 



pling in detail, we fabricated samples with a large set of SRR 
arrays (split width d = 100 nm) where d x and d y are varied in- 
dependently. We expect a blueshift with increasing density for 
all arrays since the magnetic dipoles are always transversely 
coupled. In Fig. [2 a) we plot the measured center frequency of 
the resonance versus SRR spacing for three sets of arrays. For 
square lattices d x = d y we indeed observe a continuous blue- 
shift, confirming the data for d = 80 nm in Fig.QJb). We also 
observe a blueshift when only d y is varied (d x = 500 nm). 
Remarkably, we measure a redshift when only d x decreases 
and d y is fixed at 500 nm. This result is surprising since red- 
shifts imply longitudinal coupling, which is inconsistent with 
the orientation of the magnetic dipoles. The redshift can only 
be understood by noting that SRRs also have an electric polar- 
izability in addition to a magnetic dipole liol 18 ] . The electric 
dipole moment points along the SRR base, hence allowing for 
longitudinal coupling. We implement a model that takes into 
account simultaneous electric and magnetic dipole coup ling , 
similar to the model for SRR stereodimers reported in 111 Oil . 
In this model, all magnetic dipoles couple transversely while 
electric dipoles transversely couple along y and longitudinally 
along x. We limit ourselves to electrostatic and magnetostatic 
nearest-neighbor coupling, ignoring electro-dynamic effects, 
the air-glass interface, and multipole corrections. However, 
this model captures the main physics embodied in our obser- 
vations. The coupled resonances are set by the Lagrangian 
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nm, staying below 550 nm to avoid grating diffraction in the 
range of the magnetic resonance. We measured polarization- 
resolved normal incidence transmission using the set up re- 
ported in Ref. 01311 . We illuminated a mm-sized area on the 
sample with a beam from a halogen lamp (5° opening angle) 
and used a 20 [im pinhole in an intermediate image plane to 
select the transmitted intensity from single 36 x 36 fim 2 SRR 
arrays, which we spectrally resolved by cooled Si CCD and 
InGaAs array spectrometers, and normalized to transmitted 
intensity through bare substrate. 

Fig- QIb) shows x-polarized transmission spectra measured 
on a sample with square lattices (d x = d y ) of SRRs with split 
width d = 80 nm. We observe the magnetic resonance at 
1 .4 /im only for polarization along x, as reported by [EBB, 
as well as higher order plasmon resonances at 500 nm and 800 
nm |5j|. In this Letter we focus on the magnetic resonance. 
Tracing the minimum in transmission versus SRR density in 
Fig. 01b), we find that the resonance blue-shifts as SRRs are 
brought closer. A blueshift upon increased coupling is ex- 
pected by analogy with plasmon hybridization 111, 12, 19ll . 



since the magnetic dipoles are all oriented perpendicular to the 
SRR plane, and hence transversely coupled. To study this cou- 



where L is the SRR inductance, Qi,j(Qi,j) represents the 
charge (current) on the SRR at site and where Mh and 

M e quantify the mutual inductance and the electric dipole 
coupling. Solving Eq. ([TJ for the resonance frequency at nor- 
mal incidence (feii = 0) yields 



OJ = UJQ 



\ 



1 



2«h 

"ST 



2«h 



(2) 



where ojq is the resonance frequency of a single SRR, and 
K e ,h — Me,h/L. Eq. (0 is similar to a prediction by Marques 
et al. HI for 3D SRR arrays. We fix n e = 1.04 • 10~ 21 m 3 
to match the SRR electric polarizability to the resonant ex- 
tinction cross section of 0.3 /im 2 measured by Husnik et 



al. Il20l, 21], and we set Kh = 0.67ft e . Fig. [3} a) shows the res- 
onance oj versus d x y assuming magnetic coupling only (n e — 
0). The resonance blue-shifts for decreasing d x , y in all cases 
due to transverse magnetic dipole coupling. Fig. |3b) shows 
the resonance frequency for electric coupling only (k/j = 0). 
The resonance red-shifts with increasing density unless d x is 
fixed at 500 nm. This result indicates that longitudinal cou- 
pling exceeds transverse coupling in square lattices of strictly 
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FIG. 3: (Color) Electrostatic calculation of the magnetic resonance 
frequency as a function of lattice spacing. Black curves: d x = d y . 
Blue curves: d y varies at fixed d x = 500 nm. Red curves: d x 
varies at fixed d y = 500 nm. For magnetic coupling only (a), res- 
onances always blue-shift with decreasing lattice spacing, while for 
electric coupling only (b), the behavior of the resonances for d x — d y 
changes sign with respect to (c) (all couplings). Insets in (a), (b) and 
(c) are sketches of the electric and magnetic coupling between SRRs. 
Curves in (c) are reproduced in Fig. [2] 



in-plane dipoles. This behavior is indeed observed for the 
purely electric resonance at 800 nm, at least in the regime 
m&x(d x , d y ) < 400 nm where grating anomalies [22] (asym- 
metric shoulders at 750 nm in Fig.[Tfb)) do not yet set in. 

Neither model with solely electric or solely magnetic inter- 
action is consistent with the measured shift of the 1 .4 /im res- 
onance, since we observe blue-shifts in all cases except when 
d x is varied and d y is fixed. Fig. 0c) shows the calculated 
co taking into account both electric and magnetic interactions. 
As in the data, the resonance only red-shifts when decreasing 
the distance d x at large d y . In this case longitudinal electric 
coupling exceeds the sum of transverse electric and magnetic 
coupling, leading to a net redshift. For a quantitative com- 
parison with our data we plot the shifts in Fig. |3c) together 
with the data in Fig. [2] The good quantitative agreement with- 
out any adjustable parameters confirms our interpretation that 
SRRs in metamaterial arrays show strong electric and mag- 
netic dipole-dipole interactions. These interactions are best 
quantified in rectangular arrays, since in square arrays stud- 
ied sofar ||5i|6t] partial cancellation obscures the magnetically 
induced blueshift. In our comparison we used K e /nh = 1-5 
measured in ifioll for vertically stacked SRRs. The data show 
that this ratio is also relevant for dipole-dipole coupling in the 
xy-plane, allowing a direct identification of n e and Kh with 



on-resonance electric and magnetic polarizabilities (see note 
02 IIP . It is remarkable that the magnetic polarizability ctM 
is of the same order as the electric polarizability og, as op- 
posed to the normal ordering a m < «e Mm- This conclu- 
sion is in accordance with recent estimates of Merlin M23I1 . 
that SRRs have q>m comparable in magnitude to provided 
ImeAu 3> X/£, where I is the characteristic scatterer size. 
Given the dielectric constant of Au ImeAu ~ 10, and the size 
X/£ ~ 7 of our SRRs, their LC resonances are indeed ex- 
pected to be magnetic resonances with large q>m- Into the 



visible, ImeAu rapidly decreases, causing cum to vanish M23I1 . 
as argued independently in QSQ . 

A striking feature in our transmission data (cf. Fig. [Tib)) in 
addition to the spectral shifts, is the large broadening of the 
resonance as the density of SRRs increases. In Fig. |2|b) we 
plot the measured full width at half minimum (FWHM) of the 
transmission minimum versus lattice spacing. For square lat- 
tices, the width more than doubles from 1000 to 2150 cm -1 as 
the pitch is reduced from 550 to 300 nm, while for both types 
of rectangular lattices (d x or d y fixed at 500 nm) the width 
increases from 950 to 1400 cm -1 . Such broadening was also 
noted by Rockstuhl et al. [6] for square arrays. Our extensive 
data on many rectangular and square arrays allow us to quan- 
titatively identify the source of broadening. From the outset 
it is clear that the broadening is outside the scope of Eq. £[), 
since the (Ohmic) damping rate is almost independent of cou- 
pling in any electrostatic model. Instead, electrodynamical 
radiation damping, i.e., scattering loss into the far field must 
be taken into account. As all oscillators in our sub-diffraction 
lattices are driven in phase (fcn = 0), scattered light radiated 
by all oscillators interferes destructively for all angles, except 
along the transmitted and reflected direction. Since the mag- 
netic dipoles are aligned along the incident beam, they do not 
radiate any amplitude into the fcn = directions. Hence, all 
radiation damping is solely due to the induced electric dipoles. 
For a quantitative analysis we use an electrodynamical model 
for electric poin t dip oles with a Lorentzian resonance in cke 
according to II 12[|2li 12211 . centered at 1.4 /^m and including the 
material loss rate of Au, in addition to radiation damping. We 
evaluate Eqs. (8,9) in Ref. ll22ll to predict the array transmis- 
sion. This dynamic model has no adjustable parameters, since 
the on-resonance polarizability is fixed 12111 to match the ex- 
tinction cross section of single SRRs in 12011 . We find a broad- 
ening of the collective transmission resonance that quantita- 
tively reproduces the measured broadening with decreasing 
pitch for all lattices (FWHM curves in Fig. 12b)). An im- 
portant conclusion is that the large width of the magnetic re- 
sponse commonly observed for SRR arrays fllHB] is not due 
to intrinsic loss, but is quantitatively consistent with superra- 
diant decay of the electric dipoles. The collective enhance- 
ment of the single SRR radiative linewidth, already suspected 
by |01, implies enhanced scattering and a reduction of the ab- 
sorption of the array far below the albedo of single SRRs. 

Finally we correlate the resonance broadening with the 
measured transmission T on resonance. Fig. [4] shows the 
effective extinction cross section derived from our measure- 
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FIG. 4: Effective extinction cross section per SRR derived from on- 
resonance tranmission. The black dashed line indicates the cross 
section of a single SRR (from 12011 ). The cross section per SRR is 
limited by the area d x d y of the unit cell (black and red lines). 



merits through a e g — d x d y (l — T). For uncoupled scat- 
terers we expect constant a e s equal to the extinction cross 
section er ext = 0.3 /im 2 measured for a single SRR in 1I20I1 
(dashed line in Fig. 0]), as indeed almost found in our data 
for d x = d y > 500 nm. For d < 500 nm, we measure val- 
ues for cr e ff far below er cx t indicative of strong dipole-dipole 
coupling. The collective superradiant decay (Fig.[2|b)) which 
widens the resonance reduces the extinction per element to 
remain below the unit-cell area d x d y (curves in Fig. |4j». 

In conclusion, we have measured large resonance shifts as 
a function of density in SRR arrays resonant at A = 1.4 p:m. 
These shifts are due to strong near-field electrostatic and mag- 
netostatic dipole coupling. Furthermore, we observe electro- 
dynamic superradiant damping that causes resonance broad- 
ening and an effective reduction of the extinction cross sec- 
tion per SRR. Since the data show that the response of SRR 
arrays is not simply given by the product of the density 
and polarizability of single constituents, we conclude that a 
Lorentz-Lorenz analysis to explain effective media parame- 
ters of metamaterials 'atomistically' is not valid 12411 . The fact 
that the Lorentz-Lorenz picture is invalid has important reper- 
cussions: It calls for a shift away from the paradigm that the 
highest polarizability per constituent is required to obtain the 
strongest electric or magnetic response from arrays of electric 
or magnetic scatterers. Our experiments show that increasing 
the density of highly polarizable constituents to raise the ef- 
fective medium response [5] is ineffective, since superradiant 
damping limits the achievable response. To strengthen e or //, 
we propose that one ideally finds constituents that have both a 
smaller footprint and a smaller polarizability per constituent. 
We stress that even if constituent coupling modifies e and /i, 
we do not call into question reported effective medium param- 
eters or the conceptual validity thereof per se. The effective 
medium regime only breaks down when constituent coupling 
is so strong that collective modes of differently shaped macro- 
scopic objects carved from the same SRR array have very dif- 
ferent resonance frequencies or widths. In this regime inter- 



esting physics comes into view, particularly regarding active 
devices. Specific examples are array antennas for spontaneous 
emission M1411 and Tasing spasers' B15P . where the lowest-loss 
array mode will lase most easily. 

We thank Chris Retif, Dries van Oosten, Jaime Gomez Ri- 
vas, Albert Polman and Kobus Kuipers for assistance and 
discussions. This work is part of the research program of 
the "Stichting voor Fundamenteel Onderzoek der Materie 
(FOM)," which is financially supported by the "Nederlandse 
Organisatie voor Wetenschappelijk Onderzoek (NWO)." 



* Electronic address: i.sersic@amolf.nl 
[1] V. G. Veselago, Sov.Phys. USPEKHI 10, 509-514 (1968); J. B. 

Pendry, Phys. Rev. Lett. 85, 3966 (2000). 
[2] J. B. Pendry, Physics World 14, 47 (2001); C. M. Soukoulis, S. 

Linden, and M. Wegener, Science 315, 47 (2007); V. M. Sha- 

laev, Nature Photonics 1, 41 (2007). 
[3] U. Leonhardt, Science 312,1777 (2006); J. B. Pendry, D. 

Schurig, and D. R. Smith, ibid, 1780 (2006). 
[4] D. R. Smith et al., Phys. Rev. Lett. 84, 4184 (2000); W. J. 

Padilla et al., ibid. 96, 107401 (2006); S. Linden et al., Science 

306, 1351 (2004). 
[5] C. Enkrich et al., Phys. Rev. Lett. 95, 203901 (2005). 
[6] C. Rockstuhl et al., Appl. Phys. B 84, 219 (2006). 
[7] A. N. Grigorenko et al., Nature 438, 335 (2005); V. M. Shalaev 

et al., Opt. Lett. 30, 3356 (2005); G. Dolling, M. Wegener, C. 

M. Soukoulis, and S. Linden, Opt. Lett. 32, 53 (2007). 
[8] M. W. Klein et al., Opt. Lett. 31, 1259 (2006). 
[9] E. Shamonina, V. A. Kalinin, K. H. Ringhofer, and L. Solymar, 

J. Appl. Phys 92, 6252 (2002). 
[10] N. Liu, H. Liu, S. N. Zhu, and H. Giessen, Nature Photonics 3, 

157 (2009). 

[11] M. L. Brongersma, J. W. Hartman, and H. A. Atwater, Phys. 

Rev. B 62, R16356 (2000). 
[12] A. F. Koenderink and A. Polman, Phys. Rev. B 74, 033402 

(2006). 

[13] A. F. Koenderink, R. de Waele, J. C. Prangsma, and A. Polman, 
Phys. Rev. B 76, 201403R (2007). 

[14] J. Li and N. Engheta, Phys. Rev. B 74, 115125 (2006). 

[15] N. I. Zheludev, S. L. Prosvirnin, N. Papasimakis, and V. A. Fe- 
dotov, Nature Photonics 2, 351 (2008). 

[16] S. Zhang et al., Phys. Rev. Lett. 102, 023901 (2009); E. Plum 
et al., Phys. Rev. B 79, 035407 (2009). 

[17] R. Marques, F. Martin, and M. Sorolla, "Metamaterials with 
Negative Parameters: Theory, Design, and Microwave Appli- 
cations," (Wiley, Hoboken NJ, 2008). 

[18] N. Katsarakis et al., Appl. Phys. Lett. 84, 2943 (2004).; T. 
Koschny, M. Kafesaki, E. N. Economou, and C. M. Soukoulis, 
Phys. Rev. Lett. 93, 10 (2004). 

[19] E. Prodan, C. Radloff, N. J. Halas, and P. Nordlander, Science 
302, 419 (2003). 

[20] M. Husnik et al., Nature Photonics 2, 614 (2008). 

[21] We use a cxt = Anklmct with k = 2tv/X, using the Lorentzian 
electric polarizability ao = uo K e/ (^0 — ~ iwy) implicit in 
Eq. and with radiation damping l/a — 1/ao — i |fc 3 fl2ll . 
We take 7 = 1.2 ■ 10 14 s -1 for the Ohmic damping of Au from 
P. B. Johnson and R. W. Christy, Phys. Rev. B. 6, 4370 (1972). 

[22] F. J. Garcia de Abajo, Rev. Mod. Phys. 79, 1267 (2007). 

[23] R. Merlin, Proc. Natl. Acad. Sci. U.S.A. 106, 6 (2009). 



[24] C. R. Simovski and S. A. Tretyakov, Phys. Rev. B 75, 195111 
(2007). 



